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Abstract 



We simulated a growth model in 1 + 1 dimensions in which particles are 
aggregated according to the rules of ballistic deposition with probability p or 
according to the rules of random deposition with surface relaxation (Family 
model) with probability 1 — p. For any p > 0, this system is in the Kardar- 
Parisi-Zhang (KPZ) universality class, but it presents a slow crossover from 
the Edwards- Wilkinson class (EW) for small p. From the scaling of the growth 
velocity, the parameter p is connected to the coefficient of the nonlinear term of 
the KPZ equation. A, giving A ~ p^ , with 7 = 2.1ib0.2. Our numerical results 
confirm the interface width scaling in the growth regime as W ^ X^t^ and 
the scaling of the saturation time as r ~ A~^L^, with the expected exponents 
P = 1/3 and z = 3/2 and strong corrections to scaling for small A. This 
picture is consistent with a crossover time from EW to KPZ growth in the 
form tc ^ A~^ ~ in agreement with scaling theories and renormalization 
group analysis. Some consequences of the slow crossover in this problem are 
discussed and may help investigations of more complex models. 

PACS numbers: 05.40.-a, 05.50.+q 
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I. INTRODUCTION 



Surface growth processes and deposition of thin films are of great interest due to po- 
tential technological applications (such as production of nanostructures for microelectronic 
devices, the possibility of growth of designed quantum objects, magnetic storage devices, 
among others) and due to the fundamental role these systems play in non-equilibrium sta- 
tistical physics Several models have been investigated in the last decade, most of 
them involving one kind of particle and a simple microscopic aggregation rule. The com- 
petition between different growth mechanisms have received less attention, but is essential 
to describe some practical situations, such as growth of materials designed to have specific 
electronic, mechanical or magnetic properties, which involves deposition of two or more 
chemical species. In this framework, some authors considered growth models with two kinds 
of particles and different aggregation rules Other situations involving competition 
between two growth mechanisms have also been considered 

These models usually show crossover effects from one dynamics at small times t or short 
length scales L to another dynamics at long t and large L. One typical example is Kardar- 
Parisi-Zhang [KPZ) growth at small nonlinearities The Langevin-type equation 

^ = i^V'h + ^{Vhf + r^{x,t), (1) 

known as KPZ equation, was proposed as a hydrodynamic description of kinetic surface 
roughening. Here h is the height at the position x in a rf-dimensional substrate at time t, v 
represents a surface tension, A represents the excess velocity and 77 is a Gaussian noise [jl|,|^ 
with zero mean and variance (r/ (x, t) rj (^x', t'^ ) = D6'^ (jc — x'^ 6 {t — t'). When the coeffi- 
cient A of the nonlinear term is small, a crossover is observed from linear growth (A = 0, 
known as Edwards- Wilkinson theory - EW) |]10[ to KPZ behavior. 

In discrete models, the interface width, which characterizes the roughness of the interface, 
is defined as 



W{L,t) 



(2) 



for deposition in a ci-dimensional substrate of length L [hi is the height of column i at time 
t, the bar in h denotes a spatial average and the angular brackets denote a configurational 
average). For short times it scales asW ^ and for long times, in the steady state regime, it 
saturates at Wgat ~ L°'- The dynamical exponent z = a/ (3 characterizes the crossover from 
the growth regime to the steady state regime. For systems belonging to the EW universality 
class, we have = 1/2, /?o = 1/4 and zq = 2 in d = 1 {in this paper, the subscript will 
refer to exponents of the EW theory). For systems in the KPZ class, in 0? = 1, we have 
a = 1/2, /3 = 1/3, z = 3/2 [|l],|,|10|. 

Considering the crossover from EW to KPZ scaling in ci = 1, Grossmmann, Guo and 
Grant {GGG) and Nattermann and Tang [NT) [|T2| (see also the review by Forrest and 



Toral ||T3|) proposed multiscaling relations that are equivalent to 



W{L,t) = L-f{^^,j), (3) 
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in which ~ ^1^^°- GGG also proposed that the characteristic time of crossover from EW 
to KPZ dynamics was 

tc ~ A-*, (4) 

with > 0, since the EW-KPZ crossover disappears for A = 0. Through scahng arguments, 
those authors obtained = zo/{aQ + zq — 2), which gives = 4 in d = 1. This was confirmed 
through one-loop renormalization group calculations by NT. The scaling analysis of the 
KPZ equation by Amar and Family {AF) [|14| and the assumption of Family- Vicsek scaling 
| T5| were used to show that, in the nonlinear and saturation regimes. 



W(L,t)~LVV|A|-^), (5) 

in which is a scaling function and the dependence of W on the parameters v and D of 
Eq. (|ip was omitted. A generalized scaling relation equivalent to Eq. (|^), which is a more 
general result than Eq. (j^), was also obtained by Derrida and Mallick in the context of the 
connection to the one- dimensional asymmetric exclusion model |TB|. Amar and Family |n 



have shown that the scaling form (|]) also predicts a crossover exponent = 4. 

On the other hand, all previous numerical results suggested ~ 3; for instance, GGG 
obtained this value using data collapse methods |Tl[]. Thus, it would be desirable to con- 



firm numerically the scaling properties predicted for a KPZ system in order to solve this 
controversy. 

The purpose of this work is to study a competitive growth process with EW to KPZ 
crossover, involving ballistic deposition (BD) [|T],|T7| and random deposition with surface 



relaxation (Family model) |]T8| in d = 1. In this model, incident particles aggregate to the 



deposit according to the rules of BD with probability p and according to the rules of the 
Family model with probability 1 — p. It is known that the Family model is in the EW 
universality class, while BD is in the KPZ class. This competitive model was introduced 



by Pellegrini and JuUien |T9[, whose main interest was the connection to the roughening 
transition present in higher dimensions. Although it is expected that this model is in the 
KPZ class for any p > 0, the crossover in d = 1 was not studied in detail in their original 
work and, for p ^ 0.3, effective exponents very near the EW values were obtained [|19|] . 



Here we will simulate that model in order to analyze the interface width scaling in the 
nonlinear regime, the crossover to the saturation regime and to connect the parameter p and 
the coefficient A of the KPZ equation in the corresponding continuum limit. The amplitudes 
of typical saturation times and of interface width scaling in the growth regime are consistent 
with multiscaling concepts |Tl|JI^Jl^ and refine previous numerical estimates for related 



systems. The crossover exponent = 4 follows directly from our numerical results and, 
together with the observed relation A ~ p"^'^, indicate that the crossover at small p is very 
slow. The analysis of this apparently simple problem shows that, in order to obtain reliable 
asymptotic exponents governing various quantities, it is essential to account for corrections 
to the leading terms in the scaling relations. Thus, this work may also be relevant to the 
analysis of other systems with slow crossover to KPZ scaling, whose interest increased after 
the recent debate on the problem of Fisher waves and their discrete realizations in d = 1 
dimensions pO| , pT| , p2| , p3| . For that reason, the crossover effects identified in our simulations' 



data will be discussed in detail. 
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It is also relevant to point out that a related competitive model was recently studied in 
d = 1 and d = 2 |^|, showing evidence of the asymptotic KPZ behavior. However, that 



work did not study the relation between the parameters of the discrete and the continuous 
{KPZ) model nor the scaling amplitudes that will be considered here. 

The rest of this work is organized as follows. In Sec. II, we will define precisely the 
discrete model and connect it to the KPZ equation using the scaling properties of the growth 
velocity. In Sees. Ill and IV we will present results for the interface width scaling in the 
discrete model at the nonlinear growth regime and at the steady state regime, respectively. 
In Sec. V we summarize our results and present our conclusions. 



II. THE DISCRETE MODEL AND ITS CONNECTION TO THE KPZ THEORY 

We considered a model in which particles are aggregated following the rules of BD with 
probability p or the rules of random deposition with surface relaxation (Family model) with 
probability 1 — p. In BD (Fig. la), the incident particle follows a straight trajectory 
perpendicular to the surface and sticks upon first contact with a nearest neighbor occupied 
site. It leads to the formation of a porous deposit. In the Family model (Fig. lb), the particle 
falls towards the surface along the incident column and sticks at the top of that column if 
its height is lower than or equal to the heights of the neighboring columns. Otherwise, the 
particle diffuses to the neighboring column which has the lowest height and, if two or more 
neighbors have the same height, it chooses one of them randomly. 

For p = 0, we have the Family model, which is in the EW universality class. For any 
p 7^ 0, in d = 1, we expect the BD process to change the universality class to KPZ in the 



continuum limit (see the analysis in Ref. for a related model). Then the coefficient A of 



the nonlinear term vanishes with p in the form 

A ~ p\ (6) 

with 7 > (to be estimated below). For small p and sufficiently large L, the interface width 
W{L, t,p) must scale analogously to the weak coupling regime of the KPZ theory [|13],|19|, in 
which three regimes were identified: a linear (EW) growth regime at early times {t -C tc), 
a nonlinear (KPZ) growth regime for tc t <^ t and the saturation regime for t ^ r, as 
illustrated in Fig. 2 (r is the characteristic time for the interface width saturation). 

In order to calculate the exponent 7, we considered the scaling of the interface growth 
velocity. The difference between the growth velocity in an infinitely large substrate, foo, and 
the velocity in the steady state of a finite lattice (thick films), v{L), scales as [ PB| , P7 



^v{L) = - v{L) ^ XL-''\ (7) 



with ay = 1 in (i = 1 [^. Defining 

b^{L) = Av{L) X L, (8) 



we expect that, as L — > 00, 



K{L) ^B, = BX, (9) 



4 



where S is a constant. 

In the discrete model, A varies with p, consequently by is a function of L and p which 
has a finite limiting value By{p) as L — >• oo. For very large L, Eqs. (||) and (H) shows that 
By scales with p with exponent 7. 

Simulations of the model were performed in lattice of lengths from L = 16 to L = 4096 
until the saturation regime, and in lattices with L = 2^^ = 65536 during the growth regime 
(linear and nonlinear), for several values of the probability p between p = 0.15 and p = 0.5. 
The results presented in this paper are averages typically over 10^ realizations for the smallest 
lattices [L < 256), 10^ realizations for 256 < L < 4096 and 10^ realizations for L = 65536. 
The growth velocities were calculated from numerical derivatives of the average heights of the 
deposits, with accuracies from 5 to 6 decimal places, in lattices of lengths L < 128 {L < 512 
for p = 0.15). We considered the data for L = 65536 as representative of an infinite lattice 
in the growth regime (some simulations in L = 131072 supported this assumption), and also 
obtained Voo with high accuracy. These data provided estimates of by{p,L) with accuracy 
from 0.5% to 5%. For larger lengths, poorer results were obtained due to the much smaller 
number of realizations. 

In Fig. 3 we show by{p,L) versus 1/L for the three smallest values of p considered in 
this work. The variable 1/L in the abscissa was the best choice to represent finite-size 
corrections in 6^ as L — > 00, and is related to higher order terms [l/L"^) in Eq. (0). Such 
scaling corrections have been previously observed in the analysis of small L data for BD 
and for the restricted solid-on-solid model by Krug and Meakin p6[. The corrections were 
considered in the extrapolation of the data in Fig. 3, which provided estimates of By{p) for 
several values of p (intercepts with the vertical axis in Fig. 3). 

Crossover effects may be crucial in the extrapolation procedure discussed above, and 
may severely affect the estimates of By for small p. For 0.25 < p < 0.5, four values of by 
(16 < L < 128) were well fitted by straight lines in the by{p,L) x 1/L plots (these results 
were not shown in Fig. 3, except for p = 0.25). For p = 0.2, the data for 32 < L < 256 
confirm the presence of the 1/L correction and was also used to estimate By (the estimate 
for L = 16 deviates from this trend). On the other hand, for p = 0.15, the result for 
L = 256 showed a crossover in by, which suggested calculations for L = 512. Fig. 3 shows 
that 6^,(0.15, L) slowly increases for 16 < L < 128, but decreases for 128 < L < 512. 
Consequently, the extrapolation considered only the three last points (see Fig. 3) and gave 
By ^ 0.057. However, if the extrapolation to L — >■ 00 was performed only with results for 
L < 128, then a 7% larger value of By would be obtained. Smaller values of p were not 
studied here because such crossover would appear for much larger L and, consequently, the 
extrapolations based on small systems' data would provide unreliable estimates of by{p, 00). 

In Fig. 4 we show In [-B«(p)] versus Inp using the extrapolated values of By, as discussed 
above. The linear fit in Fig. 4 gives By p^'^. Considering the error bars in By, we obtain 
an exponent 7 = 2.1 ± 0.2 (Eq. 

The large value of the exponent 7 explains the crossover effect discussed above. Since A 
decreases rapidly with p, the coefficient of the leading term in Af (Eq. 0) is small compared 
to higher order corrections (1/L^, 1/L^ etc) for small p. Thus, very large values of L are 
needed to provide reliable extrapolations with a single correction term, which prevented us 
to use values of p < 0.15 in our study. 
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III. INTERFACE WIDTH AT THE NONLINEAR GROWTH REGIME 



In the nonlinear growth regime {tc t <^ X^^L^) for sufficiently large substrates in 
d = 1, the interface width does not depend on L (weak finite-size effects). Then the scaling 
function of Eq. (j^) behaves as 

g{x) ^Cx^ , p = 1/3 , (10) 

with constant C, so that W does not depend on L, except for vanishing corrections to 
scaling. Consequently, the A-dependent scaling of W in this regime is 

W ^ C\H^. (11) 

In this section, we will verify this A-dependence through a careful analysis of simulations' 
data of our discrete model. 

However, first we will show that Eq. (|TT]) gives = 4 (Eq. ^ in a simple way, as follows. 
The crossover EW — KPZ (at t ~ tc) occurs when the scaling relation (pAf ) matches the 
EW scaling 



with C constant. Thus 



W{t,L)^C't'^\ (12) 



t * ~ Xh/, (13) 



then we obtain d) = ^„ . = , ^" . = 4. 

^ (/3-/3o) (zo-z) 

We conclude that the numerical test of Eq. (0), in particular of the dependence on 
the parameter A, may be used to test the proposal = 4. The ffist step is to extract the 
amplitude of scaling in Eq. ([TI|), which motivates the definition of the amplitude a{p,t) 
as 

a(p, t) = W{L oo, t)/t^/\ (14) 

In Fig. 5a we show a{p,t) versus for several values of p, using the data for L = 

2^^ = 65536. Different variables in the form [x > 0) were tested in the abscissa, but 
the variable 1 / 1^^^ of Fig. 5a provided the best linear fits for most values of p. The fact 
that in Fig. 5a a{p, t) is still decreasing for large t indicates the presence of a constant 
(independent of t) correction to the leading behavior in Eq. ([TTp. It proves again the 
relevance of accounting for scaling corrections in this problem, although we are not able to 
justify these corrections on theoretical grounds. 

As t ^ oo, a{p,t) converges to a finite limiting value 

A{p) = a{p, oo). (15) 

A{p) is the complete amplitude of scaling of the interface width in the nonlinear growth 
regime (Eq. |TT]). Our estimates of A{p) were obtained from linear extrapolations of a{p, t) x 
plots to t ^ oo (intercepts with the vertical axis in Fig. 5a). 
From Eq. (p^D, it is expected that the amplitude A{p) scales as A'^. From the connection 
relation (|), it is expected that 
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A{p) ~ 



(16) 



with 

5 = 7/5. (17) 

Then, the test of Eq. (|TT|) reduces to the test of Eq. (|T^ for the amphtude exponent 6. 

In order to calculate the exponent 6 in the relation (p!6D , our first step was to plot 
log [v4(p)] versus logp, but we noticed that it showed decreasing slopes as p decreased. We 
analyzed the evolution of the slopes of log [A{p)] x logp plots by calculating the following 
effective exponents for consecutive values p = p' and p = p": 

so that, as p ^ (A ^ 0), we expect that 6p 6. 

In Fig. 5b we show 6p versus p^, which gives a reasonable linear fit and indicates that 
6 = 0.7 ± 0.2. Again the variable p"^ in the abscissa is the one that provides the best linear 
fit of the central estimates of 6p, chosen among other variables in the form p^ {y > 0). In 
Fig. 5b, the effective exponents systematically decrease as p decreases, which reflects our 
previous observation of decreasing slopes in log [^(p)] x logp plots. 

Our estimates 6 = 0.7 ± 0.2 and 7 = 2.1 ± 0.2 (Sec. ^ are consistent with relation (|T7p 
with (3 = 1/3. Even considering that the error bars are large, it is relevant to notice that 
the central estimates confirm that relation exactly, which gives additional support to our 
analysis. 



IV. INTERFACE WIDTH NEAR AND AT THE STEADY STATE REGIME 

Our numerical results in the steady state regime provide additional support for the scaling 
picture proposed for the problem. 

From Eq. (^), we expect that the crossover from the nonlinear to the steady state regime 
takes place at a characteristic time r that scales as 

r ~ X-^L\ (19) 

with z = 3/2 in d = 1. In this section, our main purpose is to test the A-dependence of this 
characteristic time in our discrete model. 

The saturation time r is usually estimated using some arbitrary recipe. Here, instead of 
estimating the saturation time r (which may be defined from the time dependence of the 
interface width as it converges to the saturation value), we calculated a characteristic time 



To which is proportional to r, according to a recently proposed method That method 
provided accurate estimates of dynamic exponents for several growth models in ci = 1 and 
d = 2, including the Family and the BD models. 

First, the saturation width Wg is estimated, for fixed p and L. Then we define Tq through 

W{ro) = kWs, (20) 



7 



with fixed k {k < 1) [|^. Using the Family- Vicsek relation W{L,t) = L"f{tL~^) and 
considering that Wg ~ L'^-, we conclude that Tq ~ L^, i.e., Tq is proportional to the saturation 
time r. For the particular case of a KPZ system, Eq. (^) gives 

To ~ \'^L^'\ (21) 

Extending the procedure of previous work pS|, we considered k = 1 — 1/e = 0.6321 . . . 



in Eq. ( ]20|) to estimate tq. This value of k gave tq ~ r for BD, where r was estimated from 
the decay oiWs — W In the present model, for fixed p, we calculated the ratios tq/ L^f^ 
for several lengths L and obtained the asymptotic amplitude 

D{p) = , L ^ oo. (22) 

The extrapolation procedure follows the same lines of the calculation of A{p) from a{p,t) 
in Sec. |T|. However, only results for p > 0.2 could be obtained using data for lattice sizes 
L < 4096, since the saturation for smaller values of p is typically of EW type (r ~ L^) in 
this range of L. 

From Eqs. (|2^ ) and (^TD , we expect that -D(p) ~ A^^. Consequently, it must scale with 



p as 

Dip)r^p~^. (23) 

In Fig. 6 we show InD(p) versus Inp, with a linear fit that gives D{p) ~ p^"^-^. This result 
is consistent with the independent estimate of 7 from Eq. (H) (Sec. H). 

We also analyzed the scaling of the saturation width Ws- For lattice sizes L < 1024, 
we obtained Wg ~ L" with a = 1/2 and weak corrections to scaling. Using the data for 
L = 1024, we defined 

AWs = Wsip)-WsiO) (24) 

as the difference between the saturation width for a given probability p and the saturation 
width for the Family model {p = 0). 

In Fig. 7 we show InAiy^ versus Inp. The linear fit suggests AWg ~ p^^^, thus we 
obtain the complete form for the saturation width as 

^ (Ci + C2P^/2)^a (25) 

with a = 1/2, Ci and C2 constants. The amplitude of Ws scaling is (2^;) ^ 01) i-^-, the 
heights' fluctuations depend only on the parameters u and D of the KPZ equation (|I|), but 
not on the nonlinearity parameter A. Thus we conclude that the dependence on p in Eq. 
(p5|) is related to the dependence on p of the surface tension parameter v. when p decreases, 
the amplitude in Eq. (pSj) decreases, then the parameter v increases. Indeed, this term is 
physically expected to increase in the crossover from BD (low v) to the Family model (high 
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V. SUMMARY AND CONCLUSIONS 



We studied a competitive growth model in 1 + 1 dimensions involving two dynamics: 
ballistic deposition with probability p and random deposition with surface relaxation (Family 
model) with probability 1 — p. This model is a discrete realization of the continuum KPZ 
equation with an adjustable nonlinear coupling A related to p. At the critical probability 

= 0, the process belongs to the EW universality class, while any finite value of p drives 
the system to KPZ class. 

We established the connection between the parameters p and A as A ~ p^'^ and showed 
that W ~ p^-'^t^/^ in the growth regime. This indicates that the discrete model presents a 
very slow crossover from EW to KPZ scaling at small values of p, since the crossover time 
is tc ~ A~^ ~ This slow crossover explains the discrepancies in the effective exponents 

j3 measured in that regime in previous works [p!9|] . 

We also obtained the saturation time r ~ p~^'^L'^/^. The condition r ^ tc is necessary to 
observe the crossover to KPZ scaling, while the opposite condition leads to EW saturation 
without an intermediate KPZ growth of the interface width. A critical system size 
separates systems which present EW or KPZ saturation, and can be estimated from the 
condition r ~ tc, which gives ^c ~ A~^ ~ This large exponent proves that simulations 

in very large system sizes are necessary in order to observe all features of KPZ scaling for 
small p. 

Our results are consistent with the scaling theories for the weak coupling regime of the 
KPZ equation proposed by several authors and refine previous numerical analysis. Then 
we expect that the methods presented here may be helpful to analyze other growth models 
with slow crossovers to KPZ scaling, in which scaling theories cannot be easily developed. 
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FIGURES 



FIG. 1. (a) The aggregation rules of ballistic deposition, in which the sticking position of each 
incident particle is marked with a cross, (b) The aggregation rules of the Family model, in which 
the relaxation of incident particles to their sticking positions is indicated by arrows. The incident 
particle at the right has equal probabilities to stick at any one of the neighboring columns. 

FIG. 2. For small values of p and sufficiently large L, the interface width W{t,L,p) presents 
three regimes: a linear (EW) growth regime at early times {t -C tc), a nonlinear (KPZ) growth 
regime for tc t <^ t and the saturation regime for f^r. 

FIG. 3. by{p,L) = Av X L versus 1/L for p = 0.25 (squares), p = 0.2 (triangles) and p = 0.15 
(crosses). Solid lines are least squares fits of the data for larger L. 

FIG. 4. Log-log plot of By{p) as a function of p. The linear fit gives an exponent 7 = 2.1 it 0.2. 

FIG. 5. (a) Interface width amplitude in the nonlinear regime a{p,t) = W {L ^ oo,t) /t^/^ as 
a function of 1/i^/^ for p = 0.4, p = 0.3, p = 0.25, p = 0.2 and p = 0.15 (from top to bottom); (b) 
Effective exponent Sp versus p'^, with a linear fit that gives the exponent 5 « 0.7 as p ^ 0. 

FIG. 6. Log-log plot of the amplitude D{p) (tq/L^/"^ as L ^ 00) as a function of p. The linear 
fit gives D{p) ~ p~^ with w « 2.1 . 

FIG. 7. Log-log plot of /SWs = Ws{p) — Ws{^) as a function oi p, using data for L = 1024. 
The linear fit suggests AWg ~ p^^^. 
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